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Abstract. Planetary and lunar ephemerides are no longer used for the 
determination of inertial space. Instead, the new fundamental reference 
frame, the ICRF, is inherently less susceptible to extraneous, non-inertial 
rotations than would be a dynamical reference frame determined by the 
ephemerides. Consequently, the ephemerides are now adjusted onto the 
ICRF, and they are fit to two modern, accurate observational data types: 
ranging (radar, lunar laser, spacecraft) and VLBI (of spacecraft near 
planets). 

The uncertainties remaining in the inner planet ephemerides are on 
the order of 1 kilometer, both in relative positions between the bodies 
and in the orientation of the inner system as a whole. The predictive ca- 
pabilities of the inner planet ephemerides are limited by the uncertainties 
in the masses of many asteroids. For this reason, future improvements 
to the ephemerides must await determinations of many asteroid masses. 
Until then, it will be necessary to constantly update the ephemerides with 
a continuous supply of observational data. 


1. Introduction 

A dynamical reference frame is inherent within a given set of ephemerides, 
and the accuracy of the frame depends directly upon the accuracies of the 
ephemerides themselves. In the past, the ephemerides and their reference frame 
played a significant role in the determination of inertial space: the motions of 
the Sun, Moon, and planets were used to eliminate rotations from existing star 
catalogs. 

The situation has changed. The recently adopted IAU International Ce- 
lestial Reference Frame (ICRF) is assumed to be free from rotation; that is, 
any systematic motion of the ICRF radio sources is assumed to be negligible 
for present-day considerations. Further, modern ephemerides are automatically 
adjusted onto the ICRF through the use of ICRF-based VLBI measurements of 
spacecraft near planets. Similarly, catalogs in various wavelengths are adjusted 
onto the system of the ICRS. As such, the dynamical reference frame is no longer 
used for the determination of non-rotating reference frames. 

This paper discusses the creation of modern ephemerides and how the var- 
ious ephemeris quantities are determined by modern observational data. The 
orientation of the ephemerides onto the ICRF is shown, and the accuracies of 
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the various ephemeris quantities are discussed. Finally, it is mentioned that the 
concepts of ecliptic, obliquity, and mean equator are now obsolete in the context 
of modern ephemeris creation. 


2. Ephemerides Determine Inertial Space 

The motions of the Sun, Moon, and planets were used at one time to indicate 
extraneous rotations embedded within the stellar catalogs. This is a consequence 
of the fact that ephemeris motion, modeled with realistic equations of motion and 
fit to observational data, cannot permit a spurious rotation with respect to the 
reference system. As a consequence, the ephemerides have acted as indicators 
of systematic errors in the proper motions of star catalogs. Certainly, if the 
ephemerides were not based upon realistic equations of motion, but were the 
result of only curve-fitting instead, then the motions of the orbital planes and 
the ellipses would not be constrained by physics and the inertial basis would be 
lost. 

This may be most easily envisaged for a one-body solar system: the body 
moves in a straight line with constant velocity. Since the equations provide for 
only straight-line motion, the ephemeris of the body is automatically cast in an 
inertial frame. For a Newtonian 2-body system, the Keplerian orbit does not 
rotate in space; one simply cannot get a rotating ellipse using Newtonian, two- 
body equations of motion. For such a system, the perihelion direction acts as an 
inertial pointer; the period of an orbit, measured between successive perihelion 
passages, is determined with respect to inertial space. Similarly, for the real 
solar system, the small rotations of an elliptical orbit are due to relativistic 
effects and to gravitational perturbations from the other bodies of the system. 
However, these effects are well-known and well-modeled. 

Correct gravitational physics, along with the observational data, demands 
an inertial reference frame in order for the observations to be fit properly. If the 
observational residuals contain systematic, non-inertial errors, it is not the fault 
of the frame; the problem lies with the observations themselves or with their 
reductions. If the orbit fits the observations well, then the ephemeris is accurate 
in inertial space, mean motion and all. 

Even with fitting to only optical data, an ephemeris adjusts itself onto an 
inertial frame. The period of a planet can be determined and the motion of 
the perihelion constrained by physics; similarly, the motion of the orbital plane 
is determined by physics. Thus, a rotation in a catalog will introduce incon- 
sistencies between the physical model and the observations, albeit with limited 
accuracies. Previously, uncertainties in precession, catalog proper motions, etc., 
determined from optical observations alone, were on the order of T'/cty. 

Modern ephemerides of the inner planets, Moon, and Sun are now fit to 
highly accurate ranging observations: radar (uncertainties of 100-200 meters, 
spanning 3 decades), lunar-laser-ranging (now at the 2-3 cm level), and space- 
craft-ranging (10 meter accuracy). See Stan dish et al. (1995) and Stan dish 
(1998) for more details about the observations. Fitting to these observations 
determines all inner planet mutual distances and angles to an accuracy of about 
1-2 kilometers (1-2 milliarcseconds), and it determines the mean motions with 
respect to inertial space to the level of about 0"01/cty (equivalent to about 10 
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kilometers per century). Williams and Standish (1989) describe how well various 
orbital elements are determined by the fitting process and provide estimates of 
the uncertainties. 

Importantly, there is one part of the ephemerides that is not provided by 
fitting to ranging observations: the orientation of the whole system onto an 
external reference frame. Fortunately, this orientation is now provided by ICRF- 
based VLBI observations of the Magellan spacecraft in orbit about Venus and 
the Phobos spacecraft on its approach to Mars. These observations are discussed 
next. 


3. VLBI of Spacecraft near Planets 

VLBI measurements of a spacecraft with respect to the ICRF radio sources 
have been combined with a planet-centered spacecraft trajectory in order to 
yield a positional determination of the planet. Fitting to these observations 
automatically orients the ephemeris onto the ICRF, and since ranging observa- 
tions tightly connect all of the inner solar system ephemerides, it is the full inner 
planet system that is oriented onto the ICRF. 

A VLBI observation, taken over a relatively short time-span, is essentially a 
one- dimensional measurement. For observations taken between Goldstone and 
Madrid, where the latitudes of the two antennas are roughly the same, the deter- 
mination is almost completely in right ascension; the corresponding uncertainty 
ellipse, very long and narrow, is oriented nearly vertically in an a — 8 plot, /.e., 
the uncertainty lies along the <Caxis: a is well-determined; 8 is not. For obser- 
vations between Goldstone and Canberra, the ellipse is tilted at about 45° from 
vertical, so that the determination is a linear function between a and 8. 

The four plots in Figure 1 illustrate four VLBI observations. In each, the 
Id uncertainty ellipse is shown; the one- dimensional determination, itself, is 
the loci of points midway between the two lines of the uncertainty ellipse; the 
distance from this line to the origin represents the minimum residual of the 
possible combinations of Aa and A 8. In the upper left figure, the observation of 
Magellan at Venus was evidently taken between Goldstone and Madrid, showing 
roughly, Aci ~ — 0"003 ± 0"002. The upper right figure, between Goldstone and 
Canberra, indicates a linear combination of a and 8. The observations in the 
two lower figures are from Phobos at Mars; the wider lines of these two ellipses 
indicate the higher uncertainties for the Phobos observations. 

Table 1 lists all of the VLBI observations for Venus and Mars. The angle, 0 , 
gives the direction of the determination from the origin 90° from the uncertainty 
ellipse). The distance, is, is measured from the origin in the direction of 6; it 
gives the distance of the observation line from the origin. 

It is evident from the table that the full set of observations determines the 
orientation of the ephemerides with an uncertainty of about ± 0"001. For typical 
distances of the inner solar system, this is equivalent to an uncertainty on the 
order of 1 kilometer. 



milliarcseconds milliarcseconds 
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Figure 1. VLBI observations of spacecraft near to Venus and Mars. 
Four illustrative examples are shown; see text for explanation. 
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Table 1. VLBI observations of spacecraft near to Venus and Mars. 
The angle, 0 , gives the orientation of the observation (see the text); the 
residual, is, its uncertainty, <7, and their ratio are given. It is evident 
that the resultant orientation of the inner solar system ephemerides is 
accurate to about ±0"001. 


Venus 



0 

V ± (7 

H/cj 

1990 Sep 12 

45.6 

— 2.4± 1.0 

2.40 

1991 Dec 22 

43.2 

0.4 ± 2.0 

0.20 

1991 Dec 23 

41.9 

— 2.6± 2.9 

0.88 

1992 Mar 29 

45.0 

1.0 ± 1.8 

0.53 

1994 Apr 01 

3.1 

— 1.7± 1.8 

0.92 

1993 Nov 10 

48.8 

0.2 ± 1.1 

0.22 

1994 Feb 03 

46.0 

3.7 ± 2.0 

1.85 

1994 Mar 07 

44.4 

1.4 ± 2.2 

0.64 

1994 May 21 

-6.2 

— 2.2± 1.9 

1.15 

1994 May 22 

49.9 

7.5 ± 2.5 

2.97 

1994 May 28 

53.2 

6.6 ± 2.1 

3.06 

1994 Jun 25 

-1.0 

1.5 ± 2.3 

0.68 

1994 Jun 26 

52.3 

1.2 ± 1.9 

0.63 

1994 Jul 08 

2.7 

3.9 ± 2.6 

1.49 

1994 Jul 09 

46.0 

3.0 ± 3.2 

0.92 

1994 Jul 21 

3.4 

2.6 ± 4.0 

0.64 

1994 Aug 02 

3.5 

3.3 ± 2.4 

1.36 

1994 Aug 10 

44.4 

— 1.3± 1.3 

0.98 

1989 Feb 17 

8.2 

— 3.0± 6.2 

0.49 

1989 Mar 25 

42.3 

0.4 ± 3.3 

0.13 


Mars 
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4. Ephemeris Uncertainties 

As mentioned earlier, the major part of the uncertainties of the inner planet 
ephemerides is due directly to the uncertainties in the masses of the asteroids 
which perturb the inner planets. For the JPL ephemerides, the asteroids, 0001 
Ceres, 0002 Pallas, 0004 Vesta, 0007 Iris, and 0324 Bamberga are handled sepa- 
rately. They are each assigned a mass and their perturbations upon each of the 
major planets is included in the equations of motion of the integration program. 
Perturbation effects for 295 other asteroids are included, using masses computed 
from each asteroid’s estimated diameter and using a density according to the 
asteroid’s taxonomic class. 

The asteroid’s perturbations disturb the orbital motions of Mars and the 
Earth. Uncertainties in the asteroids’ densities and volumes lead to uncertainties 
in the masses, sometimes by as much as a factor of two. Consider covariances, 
Monte Carlo experiments using actual integrations, and sensitivity studies — all 
give approximately the same results: the uncertainties in the relative positions 
and motions of the inner planet ephemerides, introduced from the uncertainties 
of the asteroid masses, are on the order of 1 kilometer. These uncertainty esti- 
mates are also in agreement with those of Williams (1994) in his discussion of 
the asteroids’ effects on ranging to Mars. It is coincidental that the uncertainties 
from the asteroid masses are the same size as the uncertainties from orienting 
the ephemerides onto the ICRF using VLBI. 

It should be mentioned, as an aside that the uncertainties for the Moon are 
due to something else entirely. They come from tidally-induced deceleration in 
longitude which has an uncertainty of about 0"5/cty 2 . 


5. Future Improvement to the Ephemerides 

Accurate ephemerides for future dates must rely on one of two things: either the 
ephemerides must be updated continuously with current, accurate observational 
data, or there must be some significant improvement in the knowledge of aster- 
oid masses. For the latter, there are two evident methods: either the masses of 
many asteroids must be determined individually, or the method for estimating a 
mass for a given asteroid must be improved. This latter requires improvements 
in the estimation of asteroid volumes, improvements in the measuring of taxo- 
nomic classes, and improved understanding of how to assign a density given the 
taxonomic class. 


6. Conclusions 

There still is a dynamical reference frame, as in the past. It is still defined by 
the planetary and lunar ephemerides. Now, however, it is made to be coincident 
with the ICRF. The ICRF is assumed to be inertial; certainly, at least, any 
non-inertial motion of the ICRF is far less than that of the dynamical reference 
frame. As a result, the use of the dynamical frame for the determination of 
inertial space is no longer an integral part of astronomical reference frames. 

Ephemeris creation now depends upon two major types of observations: 
ranging and VLBI. For the data reductions, it is necessary to locate the observ- 
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ing stations using ICRF-based Earth orientation parameters from the IERS. 
Since all three reference frames involved in the process (ephemeris, VLBI, Earth 
orientation) are now based upon the ICRF, there is no longer an explicit use 
of the celestial equator, equinox, or ecliptic in the ephemeris creation process. 
These elements are obsolete in that context. Of course, it should be stated that 
the reduction of older observational data, still useful for the outer five planets, 
does indeed rely upon such elements. 

The accuracy of present ephemerides is limited primarily by the uncertain- 
ties in the masses of the asteroids. Future improvements to the ephemerides 
must await significant progress in the determination of asteroid masses. 
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